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Abstract. Classical reversible circuits, acting on w bits, are represented 
by permutation matrices of size 2” x 2"'. Those matrices form the group 
P(2 U '), isomorphic to the symmetric group S 2 ". The permutation group 
P(n), isomorphic to S„, contains cycles with length p, ranging from 1 to 
L(n), where L(n) is the so-called Landau function. By Lagrange inter¬ 
polation between the p matrices of the cycle, we step from a finite cyclic 
group of order p to a 1-dimensional Lie group, subgroup of the unitary 
group U(n). As U(2™) is the group of all possible quantum circuits, act¬ 
ing on w qubits, such interpolation is a natural way to step from classical 
computation to quantum computation. 


1 Introduction 

Too often conventional (every-day) computers and (futuristic) quantum com¬ 
puters are considered as two separate worlds, far from each other. Conventional 
computers act on classical bits, say ‘pure’ zeroes and ones, by means of Boolean 
logic gates, such as AND gates and DR gates [I. The operations performed by 
these gates usually are described by truth tables. Quantum computers act on 
qubits, say complex vectors, by means of quantum gates, such as ROTATOR gates 
and T gates [2]. The operations performed by these gates usually are described 
by unitary matrices. 

Because the world of classical computation and the world of quantum com¬ 
putation are based on such different science models, it is difficult to see the 
relationship (both analogies and differences) between these two computation 
paradigms. In order to remedy this problem, in recent years, NEGATOR gates and 
PHAS0R gates have been developed as basic building-blocks of quantum circuits 
III- In the present paper, we bridge the gap between the two sciences by 
deriving the NEGATOR gates in an alternative way. The common tool we have 
chosen for describing both reversible computation and quantum computation is 
the matrix representation. Classical reversible circuits iTj, acting on w bits, are 
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represented by permutation matrices of size 2 W x 2 W , whereas quantum circuits 
[2], acting on w qubits, are represented by unitary matrices of size 2 W x 2 W . 
Invertible matrices form a group under the operation of matrix multiplication. 
The matrix group consisting of permutation matrices is a subgroup of the group 
of unitary matrices. In the present paper, we show how to enlarge the subgroup 
to its supergroup, in other words: how a classical computer can be upgraded to 
a quantum computer. It is surprising that a mathematical tool for this tour-de¬ 
force is the good-old polynomial interpolation formula of Lagrange. By Lagrange 
interpolation between two (or more) permutation matrices we indeed obtain an 
infinity of unitary matrices. 

2 The Lagrange interpolation 

We consider the n x n unitary matrix q. The only restriction is the finiteness 
of its order. In other words: we assume that q p equals the n x n unit ma¬ 
trix u and that none of the matrices q s with 0 < s < p equals u. Thus the 
set {q,q 2 ,q 3 ,...,q p ~ 1 ,q p } (with q p = q° = u) constitutes a finite matrix group 
isomorphic to the cyclic group Z p of order p. 

We construct a 1-dimensional Lie group which constitutes a smooth inter¬ 
polation between the p matrices q 3 . For this purpose, we are inspired by the 
interpolation formula of Lagrange: if, for the values x = Xj, a function y(x) 
evaluates to yj, then the polynomial 



automatically satisfies p(xj) = y j, for all Xj. The function p(x ) connects the 
points ( Xj,yj ) in an analytic and smooth way. By analogy, we construct the 
matrix 



( 1 ) 


where the constant w is the primitive p th root of unity, i.e. e* 27r / p . Here, indice^l 


(j and k ) run from 0 to p — 1. The reader can easily verify that indeed we auto¬ 
matically have m(j 2tt/p) = q 3 . In particular we have to(0) = q° = u. Expression 
m constitutes a generalization of the case p = 2, discussed in Appendix A of 
[5]. Indeed, for p = 2, eqn dTJ) recovers the expression found in [5j: 



( 2 ) 


Rewriting m as 



(3) 


o 


1 Also in the rest of the paper, each time the limits of a summing index are not 

specified, they equal 0 and p — 1. 
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the coefficients nij(6) (known as the Lagrange basis polynomials or Lagrange 
fundamental polynomials [6] ) have the property that m,j(k 2ir/p) equals 1 for 
k = j but equals 0 otherwise. This is illustrated for p = 3 in Figure [1] where we 
see the modulus squared of the three coefficients 


m 0 {d) = 
mi(0) = 
m 2 {0) = 


1 — CO 

e ie - 1 

CO — 1 

e ie - 1 
co 1 — l 


1 

- u 2 


, ,2 

e 

— UJ 

<jj ■ 

-CO 2 

if) 


e 

— UJ 

CO 2 

• 

— UJ 


where co is the primitive cubic root of unity, i.e. e i2,r / 3 = —1/2 + * a/3/2. 



Fig. 1. The moduli squared of the three coefficients in rrij(9)q : ’, the La¬ 

grange interpolation between the three unitary matrices q°, q 1 , and q 2 . 


There exist various different expressions equivalent to 0: 

m(0) = - ^ / I]V e “ 

^ i 

m(0) = - (e ipe - 1 ) £ 


e l6 — oji 


i(0) = - E^E 


co~ rj e irB 


(4) 

(5) 

( 6 ) 


Appendix A describes how they are deduced from ©. Eqn (|4]) is more compact 
than eqn m, as it contains only one instead of two products. Expression m 
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corresponds with the so-called ‘first barycentric form’ of the Lagrange interpo¬ 
lation formula. It has the advantage that it contains no product at all. However, 
it has a disadvantage with respect to 0: in eqn 0 we clearly see that to is a 
polynomial (of degree p — 1) in the variable e l6 \ in eqn (0 this fact is hidden. 
The expression 0 constitutes a good compromise between 0 and 0: we still 
recognize the polynomial nature of the expression, while the sum of products in 
0 is simplified to a sum of sums in 0. 

In Appendix B, we prove that m{9) has the property 


m{6i)m{9-2) = m{9i + 9 2 ) , 


such that the multiplication of two m{9) matrices yields a third m(9) matrix. In 
particular we have m{9)m{—9) = to(0) = u, such that m(9)~ 1 = m(—9). Hence, 
all conditions are fulfilled to say that the set m{9 ), together with the operation 
of ‘ordinary matrix multiplication’, forms a group. It is a continuous group: a 1- 
dimensional Lie group. In Appendix C we demonstrate that the matrices m{9) 
are unitary. We thus may conclude that the group is isomorphic to the unitary 
group U(l). 

Any 1-dimensional Lie group m(9) has a generator: 


1 dm 



With the help of 0 we find: 




With the help of 0 we find an alternative expression for this n x n matrix: 



For the case p = 2, both expressions simplify to g = (q° — q 1 )/ 2, a result that 
can also be retrieved directly from 0. 

3 Examples 

The n x n unitary matrices form the unitary group U(n). Within this infinite 
group figures the finite subgroup P(n) of all n x n permutation matrices. We 
choose as matrix q of the previous section, one of these n\ permutation matrices. 
Such matrix generates a finite set of matrices {q, g 2 ,..., q p }. This set is a subgroup 
of P (n), isomorphic to the cyclic group Z p of order p. The value of p depends 
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on the particular choice of the permutation matrix q. The minimum value is 1 
(for the trivial choice q = u); the maximum value is L(n), where L denotes the 
Landau function [7j. 

The group P(n) is isomorphic to the symmetric group S„. The cycle graph 
of a finite group depicts its cyclic subgroups. By convention, only the primitive 
or maximal cycles (i.e. those cycles that are not subsets of another cycle) are 
shown. 

As a first example, we investigate P(2), i.e. the group of all 2 x 2 permutation 
matrices. It contains two elements, forming a single 2-cyclc: 


• the two matrices q 




= ?' 3 


the matrix q representing the NOT gate. After @, the interpolation between q° 
and q l is 


m{8) 


1 (l + e ie l-e i0 \ 

2 Vl-e ie l + e ie ) ’ 


called the NEGATOR gate [8] . It is a quantum gate, generalization of the classical 
gates to( 0) (i.e. the IDENTITY gate) and m( n) (i.e. the NOT gate). We note that, 
for any value of the angle 9 1 both row sums and both column sums of m{9) are 
equal to 1. The generator of m(9) is 



a matrix with all line sums equal to 0. 

As a second example, we take P(3), i.e. the group of all 3 x 3 permutation 
matrices. It consists of 3! = 6 elements, ordered in four maximal cycles: 


three 2-cycles: 

• the two matrices q 

• the two matrices q 

• the two matrices q = 



and q 2 = 


and q 2 = 


and q 2 = 



— one 3-cycle, consisting of the three matrices 


/ 010 \ 

9 = 001 , q 2 

\ioo ) 


/ 001 \ 

110 0, and q 3 

\oio ) 


100 ' 

0 10 I = q° 
0 0 1 
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Fig. 2. Cycle graph of the group P(3) of 3 x 3 permutation matrices. 


There exist no longer cycles, as L(3) = 3. Figure [2] displays the graph. 

We now examine the 3-cycle in detail. After and (O, respectively, its 
unitary interpolation is 


m(9) = - (e i3B - 1 ) 


e lf) - 1 


T + 


e* 9 — to 


e' e - w 2 q 


= \ [ ( e l2e + e ie + 1 ) q° + {uje i2S + w 2 e* 6 + 1) q 1 + ( oj 2 e l29 + we“ + 1) q 2 } 

O 


1 


x 2 + X + 1 UJX 2 + 0J 2 X + 1 UJ 2 X 2 + OJX + 1 


= - I UJ~X + UJX + 1 X + X + 1 UJX + UJ X + 1 
, UJX 2 + UJ 2 X + 1 UJ 2 X 2 + UJX- 1-1 x 2 + X +1 


where a; is a short-hand notation for e lB and uj again is the primitive cubic root 
of unity. As expected, we have m{ 0) = q°, m{ 2tt/ 3) = q 1 , m( 47r/3) = q 2 , and 
m(27r) = q° again. 

We notice that all six line sums (i.e. three row sums and three column sums) 
of m(9) are equal to 1. Therefore, we say that m belongs to the subgroup XU(3) 
of U(3), described in detail in Appendix B of !8). Here, XU(rt) is the subgroup 
of U(n) consisting of all n x n unitary matrices with all 2 n line sums (i.e. n row 
sums and n column sums) equal to unity m m 0- This is no coincidence. In 
Appendix D, we prove that, if q is a unitary matrix from XU(n), then also the 
interpolation is a member of XU(n). 
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For the generators of the four cycles we find: 


1 

2 


/0 0 0 \ 


\ 0 0 0 / 


1 

2 



^ / 3 w — 1 w 2 — 1 \ 

and - la ; 2 — 1 3 a; — 1 ] 

ya; — la ; 2 — 1 3 ) 


respectively. Each has all line sums equal to 0. Separately, each of these gener¬ 
ators generates a 1-dimensional subgroup of XU(3); together they generate the 
full 4-dimensional group XU(3), subgroup of the 9-dimensional group U(3). 


4 Classical and quantum computing 


Because classical reversible circuits [T] acting on w bits, are represented by ma¬ 
trices from P(2 W ) and quantum circuits [2] acting on w qubits, are represented 
by matrices from U(2 1 "), we have special attention for the case n = 2 W . 

As an example, we investigate the group P(4), representing all possible re¬ 
versible circuits with 2 bits at the input and 2 bits at the output. There exist 
4! = 24 such circuits. Because P(4) is isomorphic to S 4 , we investigate the cycle 
graph of this symmetric group. It consists of six 2-cycles, four 3-cycles, and three 
4-cycles. No longer cycles exist, as L( 4) = 4. We have a closer look at one of the 
4-cycles: 



/o 1 0 o\ 
0010 

2 

/o 0 1 0 \ 
0001 

3 

/o 0 0 1 \ 
1000 

4 

/I 000 \ 
0 100 

q = 

0001 

, q = 

1000 

, q = 

0100 

, q = 

0 0 10 


^1000^ 


^010 oy 


^0 0 1 o^i 


^ 0001 / 


For the interpolation between u , q , q 2 , and q 3 , we find 


m(6) = 


1 

4 


/1 + x + x 2 + x 3 1 — ix — x 2 + ix 3 1 — x + x 2 — x 3 1 + ix — x 2 — ix 3 

1 + ix — x 2 — ix 3 1+ x + x 2 + x 3 1 — ix — x 2 + ix 3 1 — x + x 2 — x 3 

1 — x + x 2 — x 3 1 + ix — x 2 — ix 3 1 + x + x 2 + x 3 1 — ix — x 2 + ix 3 

y 1 — ix — x 2 + ix 3 1 — x + x 2 — x 3 1 + ix — x 2 — ix 3 1 + x + x 2 + x 3 


We note that the matrix Q = q 2 not only belongs to this 4-cycle, but also is 
member of a 2 -cycle (although not a maximal 2 -cycle): 


Q = 


fooio\ 
0 0 0 1 
10 0 0 
\010 0 / 


Q = 


(l 0 0 0 \ 
0 10 0 
0 0 10 
\0 0 0 1 / 


= Q° 


For the interpolation between u and Q , we find the Lagrange interpolation 


M{B) 


(l + x 


1 

2 


0 

1 — x 


\ 0 


0 

1 + x 

0 

1 — x 


1 — x 0 \ 

0 1 — x 

1+ x 0 

0 1 + x ) 
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/0 1 0 0 \ 
0 0 10 
0 0 0 1 
\1 0 0 0 / 



/I 0 0 0\ 
0 10 0 
0 0 10 
\0 0 0 1 / 


/0 0 0 1 \ 
10 0 0 
0 10 0 
\0 010/ 


Fig. 3. One maximal 4-cycle and one non-maximal 2-cycle of the group P(4) of 4 x 4 
permutation matrices. 


different from m(9). Figure [3] displays both the 4-cycle and the 2-cycle. 

Although we have M( 0) = m(0) = u and M(tt) = m(ir) = Q 1 we have 
M(tt/2) ^ 77t(7t/ 2). Indeed, m(7r/2) = q is a permutation matrix, whereas 
M(tt/2) is a complex unitary matrix: 


m( 7t/2) 


/0 1 0 0\ 
0 0 10 
0 0 0 1 
\1 0 0 0/ 


and M(tt/2) 



fl+i 

0 

1 — i 

0 \ 

1 

0 

1 + i 

0 

1 — i 

2 

1 - i 

0 

1 + i 

0 


\ o 

1 — i 

0 

1 + */ 


They represent the circuits 


-- 0 - 

—0—•— and - > 

repectively. The former is the cascade of a controlled NOT and a NOT; the latter 
is a square root of NOT (a.k.a. a V gate). The former is a classical circuit; the 
latter is a quantum circuit. Both circuits are square roots of the classical circuit 


with matrix representation Q. We close this section by remarking that the above 
two matrix sets m(6) and M(6), both 1-dimensional interpolations between the 
unit matrix and the matrix Q , illustrate that a unitary interpolation between 
two unitary matrices is not necessarily unique. 
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Efficient circuit design applying V gates, controlled V gates, W gates, or con¬ 
trolled W gates !.10j m, can be interpreted as using cycles of dyadic 2x2 unitary 
matrices [12l as matrix q of Section 2, rather than permutation matrices. 

5 Conclusion 

By Lagrange interpolation, it is possible to embed a finite cyclic group in a 1- 
dimensional cyclic Lie group. By applying this technique to a cycle of an n x n 
permutation matrix, one obtains a 1-dinrensional subgroup of the unitary group 
U(n). In this way, we can bridge the gap between classical reversible computation 
(represented by permutation matrices) and quantum computation (represented 
by unitary matrices). 
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A Alternative expressions for the Lagrange interpolation 

Assume the numbers w°, w 1 , w 2 , and w p_1 are the p solutions of the equation 
x p —1 = 0. Hence: 

x p — 1 = (x — uj k ) , 
k 


such that 


n (*-“>*) 

k^j 


X p -1 

x — oV 


(7) 


We apply this result twice: 


— With x = oji e le and subsequently e -» 0 we obtain 




uj k ) = lim 
6—^0 


uj p ^e ipe - 1 
— uP 


lim 

e—>0 


e ipe - 1 
w-J (e i£ - 1) 


1 ipe p 

uji e—>o ie tjji ’ 


such that m becomes ©• 
— With x = e lS we obtain 


\[{e i6 -u: k ) 

k^j 


e ipd - 1 

e lS — ojj ’ 


such that © becomes ©. 

Eqn (0) is obtained from eqn (0) by computing w- 5 II — ujk )- We find its 
value as follows: 


^ u~ rj x r = J2 ^~°x) r 

r r 

_ (u- j x) p - 1 
uj-ix — 1 
_ uj-^xP - 1 
w~ix — 1 
x p - 1 
uj-ix — 1 


= oj j JJ(a; — uj k ) . 
k^j 
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B Proof of the group structure of the Lagrange 
interpolation 

Assume a general operator (e.g. in matrix form) a of the following form 

p -1 

S = £ a iQ 3 ’ 
i=o 

with a; € C elements over the complex field, and q p = 1. From now on, we 
will omit the operator hat notation on the q operators. First, a straightforward 
rearrangement shows that 

Lemma 1 Given two operators a and b, then the product operator 

c = ab 


has coefficients 


p -1 


— ^ ^ ambi—m. d“ ^ ^ . 


m —0 


m=z+l 


Note that, in Cj, the indices of the coefficients of a and b add up to i in the first 
summation and to p + i in the second summation. We take the summations to 
be strictly increasing, such that the second term of c p _i vanishes by definition. 

Proof. The proof is based on a simple rearrangement of terms. Explicitly, one 
can write 

p— 1 p—1 

«* = ££ affi k q j+k , ( 8 ) 

.7=0 k —0 

and break this double summation down in three different regions, corresponding 
to respectively regions (A), (B) and (C) in Figure U 


p—2p—j—2 p —1 p —1 p—1 

ab = ^2 £ ajbkq j+k + ^2 a jbp-j-i<f~ 1 + ^2 a jbkq j+k ■ 

j—0 k— 0 j=0 j— 1 k—p—j 

A change of dummy summation indices {i = j + k,l = ^(k — j)} givet@ 


P-2 


P _1 
2 2 




2p-2 P-1- 

+ £ £ 91 ■ 


i=P l=±-p+l 


Note that l can take half-integer values. 


2 
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♦ . 

iJ 



Fig. 4. Decomposition of the summation in eq. lf8)l . 


The index i is always larger (or equal) to p in the last summation, such that 


P-2 


p_i 
2 2 


a i-i b i+i qi + X a i -h-ih-h+i * 


p- 1 


i=° 2=-* 
p- 2 f-i-’ 


+ X X qn 


m=0 2=i£-£+l 


where q p = 1 has been used. Reintroducing dummy summation index j = | — l 
in the first summation, while defining it as j = | — i — l in the second, and 
j = I + T — l in bist, we obtain 

p—2 2 p—1 p—2 p—1 

= X X a j b i-j < f + X + X X a j b P+i-J<l l ■ 

i= 0 j—0 j=0 2=0 j=2+l 

This leads us to 


p-i 


aS = X 


2=0 


X a ™ 6 2—r 


P-1 

^ ^ Um^p+2—r; 


.m=0 


271 = 2+1 


which completes the proof of the lemma. 

The 1-parameter group structure of the Lagrange interpolation operators can 
now be proven. As demonstrated in Appendix A, there are multiple represent a- 
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tions for the rrij coefficients in (|3)l . We will employ here expression d5j) : 


i(0) = ~ 


1 e ipS - 1 


p e™ — uj 3 

Lemma 2 Given two Lagrange interpolation operators fh( 6 1 ) and 101 ( 62 ) with 
coefficients 

1 e ip8k - 1 


irij(9k) = 


to 3 for k € {1, 2} , 


p e ™k _ 

then the product operator a = 01 ( 61 ) 171 ( 62 ) has coefficients 

l e ip(6 1 +9 2 ) _ 1 


3 p e i(9i+92) _ 

Proof. The previous lemma states that 


UJ J 


3 1 e ipei - 1 1 e ipB2 - 1 

u = E ~ w - —-— u 3 


m =0 


p e lBl — u) m 


p e 1 ® 2 — uj 3 


—m 


1 e ipBl - 1 1 e ip62 - 1 

+ V - - u> m - 

z -—' p e lBl — uj m p 

m=j +1 y 


Using u; p = 1, this can be shortened to 


1 


= -=■ (e lpBl - l)(e lp ” 2 - 1W 

p 2 


P-1 


1 


V 

' pt'Q \ 


,p+j-m 


1 


m —0 


jm g id 2 _ ffi-r 


The summation in © can be reduced using partial fractions: 
p - 1 1 1 

E 

m =0 


s i6»i _ e id 2 _ w j-r 


1 P - 1 1 

L v_- 

iO 2 Z—/ _ 


0T 


e by z ' e 

m—0 


w m w m _ w j e -*02 


1 




ur 


gi(6i+e 2 ) _ kp 


' u> m - UJ 3 e~ lti2 


—i9o 


1 


e i(0i+e 2 ) _ 


p- 1 r 

E 

m—0 

P- 1 Tgi^! _ gift + ^ w m _ ^jg-ift + w i e -^2 

77T m I - 


E 

m—0 L 


1 


e i(ft+ft) _ fjjj 


p- 1 1 

s iei y 

Z_—/ pWp _ 


u m -uie- i9 * 
p -1 


— u^e 202 


a; r 


E 


1 


m—0 m=0 

This formula can be even further reduced by realizing that 

e n<*-^>=x n<*-"‘>=£(*'-1)= 


oj 3 e~ ie2 _ 


m—0 k^m 


px 


.P-1 


(9) 


( 10 ) 


( 11 ) 
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and thus 

p -1 


1 ^ ) 


y = y 

' x — U> m ' 


1 ^ TT, t, 


ru(*- wfc ) iifc(*-w fc ) 


e n <*-+= s+r • 


Inserting this relation in both sums of G0. we gather 

p-i 


1 


V _- 


1 


m—0 


e *e 2 _ 


1 


a iOi P e 


i(p-l)8i . T¥j ,(p-I)i e -*(P- I ) e 2 


e i(8 1 +e 2 ) _ W gipSi _ 2 

p 1 _ gip(01+02) 




U} pj e -ip8 2 _ X 


gi(0i+0 2 ) _ w j (gipSi _ X)(l _ gjpS2) ’ 

in which we have used that oj p = 1. Inserting this result in eq. © leads to 

X gip(0i+0 2 ) _ x 

a, =- a , -- w 7 = m,-(0i + 0 2 ) ■ 

3 p e *(ei+e 2 )_ w j ■ n 

This means that m(0i)m(02) equals to(0i + 02), which is what had to be proven. 

C Proof of the unitarity of the Lagrange interpolation 

Because of Appendix B we have the matrix equality m(0i)m(02) = m(01 + 02), 
and hence, in particular, m(0)m(—0) = u or m(0) _1 = m(—6). We compute 
m(—0) from © by changing from indices j and fc to indices J = p — j and 
A' = p — k: 

r m V- n fc ^-(e _< ®-w fc ) 

™ -0 = > o” , . -*7 


rwe-* e -u+-*) 


E llis+jl 

= TT-7^=7-~ <7 


V—J 


= £ 


7 TIayj( w J -w ■ ftr ) 


7 rwo/-^) 


(9 f )‘ 


= [m(0)]t, 


where we took into account that q is unitary: q 1 = gb 

We can conclude that m(0) _1 equals [m(0) ]t and thus that m(0) is unitary. 
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D Proof of the XU property of the Lagrange interpolation 

For the classical Lagrange interpolation, the sum of the Lagrange fundamental 
polynomials is equal to unity, a property known as the first Cauchy relation [B| 
|9j. We prove here that a similar property holds for the matrix interpolation 

CD-©: 

Y^rrijid) = 1 . 
j 

Indeed, by applying eqn (j4j) , we find: 

= - IK— fc ) 

3 3 k^j 

= ~ ~ I - wj ) “ x ] JJ(a: - w fc ) 

^ 3 k^j 

=-1( x - wi ) ik* - ^\ y, m*- wfc) • 

3 k^j j kjij 

Applying © for the former sum and m for the latter sum, we obtain 

- ( x p - 1) + x p 


and thus unity. 

If a matrix q belongs to the group XU(n), then all q 3 belong to XU(n). 
Because the row sum of a sum of matrices equals the sum of the row sums of 
the constituent matrices, the row sum of m = nijq 3 equals times 

the unit row sum of q 3 and thus equals JA 171 j ■ Because of the above Cauchy 
relation, this quantity is equal to 1. Hence, the matrix m has all n row sums 
equal to unity. We obtain a similar result for the n column sums. Hence, the 
interpolation matrix m belongs to XU(n). In fact, the matrices m(6) constitute 
a 1-dimensional subgroup of the (n — l)-dimensional group XU(n). 


